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Abstract. We study parameter estimation of supermassive black holes in the 
range 10 s — 10 s Mq by LISA using the inspiral full post-Newtonian gravitational 
waveforms, and we compare the results with those arising from the commonly 
used restricted post-Newtonian approximation. The analysis shows that for 
observations of the last year before merger, the inclusion of the higher harmonics 
clearly improves the parameter estimation. We pay special attention to the source 
location errors and we study the improvement on the percentage of sources for 
which we could potentially identify electromagnetic counterparts. We also show 
how the additional harmonics can help to mitigate the impact of losing laser links 
during the mission. 



PACS numbers: 04.25.Nx, 04.80.Nn, 95.55.Ym, 97.60.Lf 

1. Introduction 

LISA will be an astronomical observatory of unprecedented versatility and range [flU]- 
Among the wide range of different LISA sources, the observation of supermassive black 
holes, in the range 10 5 — 10 8 M Q , merging in galaxies at all distances will address many 
of LISA's science objectives. In particular, they will provide valuable information 
about the mechanism of their formation [21 2] and they will also serve as laboratories 
for fundamental tests of gravitational theory H>1 [51 [El 03 HH]- In addition, since 
many of these supermassive black hole mergers are likely to have electromagnetic 
counterparts, it will be possible to constrain the values of cosmological parameters 
by combining the gravitational wave and electromagnetic observations [H [HI [12J, [13] . 
Moreover, using the distance-redshift relation from many supermassive black holes, 
LISA will be able to put interesting constraints on the equation of state of dark 
energy [14j . The real impact of these observations will depend on how accurately the 
source parameters can be estimated. 

Supermassive black hole binaries are long lived sources in the LISA band. The 
whole coalescence of the compact binary system is usually divided into three phases: 
the adiabatic inspiral, the merger and the ringdown. Most of the signal-to-noise ratio 
(SNR) accumulates during the last days of the coalescence, but one critically relies on 
long integration times to disentangle the source parameters, in particular to resolve 
the source position in the sky and measure its luminosity distance |15j . 
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In this paper we study parameter estimation of supermassive black hole binary 
systems in the final stage of inspiral (ignoring the merger phase and the ringdown) 
using the full post-Newtonian gravitational waveforms and we compare the results 
with those arising from the commonly used restricted post-Newtonian approximation, 
complementing the results presented in [16l [17l [18] 

The rest of the paper is organized as follows. In section [2] we describe the 
waveform model we employ and spell out the assumptions on which our analysis is 
based. In section [3] we present the results exploring the vast parameter space, paying 
special attention to the source location errors, and we also compare the results in 
the hypothetical case in which LISA could lose two laser links and would operate as 
a single interferometer. Section [4] concludes with a summary and pointers to future 
work. 



2. Observation of supermassive black hole inspiral systems by LISA 

The coalescence of binary black holes is commonly divided into three successive 
epochs in the time domain: inspiral, merger and ringdown. During the inspiral the 
distance between the black holes diminishes and the orbital frequency sweeps up. 
The waveforms are well modeled using the post-Newtonian approximation to general 
relativity (see [19] and references therein). Eventually the post-Newtonian description 
of the orbit breaks down, and the black holes cannot be treated as point particles any 
more. What is more, it is expected that they will reach the innermost stable circular 
orbit (ISCO), at which the gradual inspiral ends and the black holes begin to plunge 
together to form a single black hole. This is referred as the merger phase. At the end, 
the final black hole will gradually settle down into a Kerr black hole. 

The inspiral post-Newtonian (PN) waveforms in the two polarizations h+ and 
h x , take the general form 

W 2 H^ + v 2 H^ x + v^H^ + ...}, (1) 

where we have set G = c = 1, v = (A/w) 2 ' 3 is the PN expansion factor, oj is the orbital 
frequency, Dl is the luminosity distance to the source, and M and rj are the observed 
total mass and the symmetric mass ratio respectively. The explicit expressions for 
^m/2) can k e f ounc j j n [5T] and they include contributions from several harmonics 
of the binary's orbital motion. Equation (fT]) corresponds to the so-called full waveform 
(FWF). If one neglects all amplitude terms except the leading Newtonian quadrupolc 
one, but keeping the phase to some specific PN order, i.e., keeping only ffi? x and 

throwing out the rest H^ 71 ^ 2 ' for m > 0, it becomes the restricted waveform (RWF). 

For the supermassive black hole binary inspirals most of the SNR accumulates 
at frequencies / < 10 mHz, so it is adequate to use the low- frequency approximation 
to the LISA response function derived by Cutler |15j . In this approximation, LISA 
can be regarded as two independent Michelson interferometers, and the strain h{t) 
produced by a gravitational wave signal becomes 

fe W( t ) = ^ \F®(t)h+(t) +F®(t)h x (tj\ , (2) 



2 

c • 

labels the two independent Michelson outputs 



(i) (i) 

where _FJ_ and F x are the time-dependent antenna pattern functions and the i =1,11 
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The total noise that affects any LISA observation is given by the superposition of 
instrumental sources, S!^ st (f), and astrophysical foregrounds of unresolved radiation, 
Sn" ! (f)- I n this paper we use the same analytical expressions given in [53] 

ST(f) = 6-12 x 1(T 51 r 4 + 1.06 x 1(T 40 

+ 6.12 x 1CT 37 f 2 Hz -1 . (3) 

For the confusion noise we consider only noise from short-period galactic and 
extragalactic binaries (due to white dwarfs binaries), assuming they are all 
unresolvablc, and we ignore the effects of captures of compact objects. 

In this paper we employ the Fisher matrix approach to study the problem of 
parameter estimation for supermassive black hole inspirals, paying especial attention 
to the improvement in errors estimation using the FWF in comparison to the RWF. 
The analysis performed is based on the following assumptions: 

• We consider the last year of the inspiral phase of the coalescence. We terminate 
the signal when the binary members are separated by a distance 6M and we also 
impose a low-frequency cut-off to the instrument at 5 x 10 -5 Hz. 

• For the waveform model, we restrict ourselves to circular orbits and we take care of 
spins contributions only in the waveform phase, ignoring spin-induced precession 
of the orbital plane. Moreover, we focus our attention to the cases in which spins 
are negligible. 

• We approximate the waveform at the 2PN order, both in amplitude and phase, 
considering up to six harmonics in the case of FWF. These waveforms depend on 
eleven independent parameters: 

A = {cos 6 N ,(f> N , cos 6 L ,<p L An D L ,t c ,<j) c , ft, a,]nMMlA , (4) 
four angles defining the source position and orientation, the luminosity distance, 
t c and <pc ar< 5 the time and phase at coalescence, (3 and a are the spin-orbit and 
spin-spin parameters, and two mass parameters. 

• We consider sources at rcdshift z = 1 in a flat Universe described by the 
cosmological parameters Ho = 71 km s _1 Mpc" 1 , Q m = 0.27 and fl\ = 0.73. 

3. Results 

The results we present here complement and extend those presented in [THl [T7] [TS] . 
These papers studied the effects of adding higher post-Newtonian order corrections to 
the amplitude, analyzing how the SNR and the measurement errors evolve with the 
total mass of the system for different values of the mass ratio. Their conclusions were 
the following: 

• For the SNR, the contribution of the second harmonic dominates for systems 
with M < 4 x 10 7 Mq, but for higher mass systems the second harmonic is no 
longer visible and then the contribution of the higher harmonics becomes relevant, 
increasing LISA's mass reach for supermassive black holes. 

• For the measurement errors, the use of higher harmonic terms improves the 
parameter estimation, not only for the most massive systems, but for any binary 
with total mass higher than 10 5 M Q . The angular resolution has gain factors 
from 2 — 3 up to one order of magnitude when they go to higher masses, and the 
estimation of masses has even better improvements, with gain factors of at least 
one order of magnitude. 
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Table 1. 10%, 50% and 90% levels of the cumulative probability distributions of 
SNR and measurement errors for different pair of masses using FWF and RWF, 
together with the gain factors computed as {5A r iJ 5 Q%) FWF /(SA r i? 5 Q%) RWF and 
10(( i 50%>rwf-(i 5 o%)fwf) i for the SNR and measurement errors, respectively. 
The probability distributions are obtained by considering LISA as a combination 
of two independent Michelson interferometers. 
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In this paper we analyze in more detail the effects of the higher harmonics 
corrections in the distributions of the parameter estimation errors, interpreting the 
results in a complete new way that provides more useful information about the impact 
of using the higher harmonics in modeling supermassive black holes. In particular we 
pay special attention to the source location errors and we study the improvement 
on the percentage of sources for which we could potentially identify electromagnetic 
counterparts. We also study how the additional harmonics in the FWF can help to 
mitigate the impact of losing laser links during the mission, therefore extending the 
previous results. 

For this analysis, we perform extensive Monte-Carlo simulations in order to 
extract general conclusions of which is the real impact of using these higher harmonics 
corrections, since the results can vary significantly from source to source. As we did 
in a previous paper [18] . for different pair of masses, we consider an ensemble of 1000 
fiducial sources all at rcdshift z = 1 (which, according to our Universe model, sets 
the luminosity distance to Dl = 6.64 Gpc), with zero spins (3 = a = 0, and we 
select randomly the four geometrical angles (On, 4>n, 6l and 4>l) from an uniform 
distribution in cosOn, 4>Ni cos9l and <Pl- 

The probability distributions of SNR and measurements errors for observations 
of the final year of supermassive black hole binaries, by considering LISA as a two 
independent Michelson interferometers, can be found in figure 11 of [IS]. Since LISA's 
measurement errors span several orders of magnitude with just changing the sky 
location and orientation of the source (see e.g. figures 7-10 of [18]), it is very important 
to characterize these distributions properly. For this reason, we provide in table [T] the 
10 %, 50 % and 90 % levels of cumulative probability distributions, for the SNR and the 
measurement errors, for four particular pair of masses using the RWF and the FWF, 
together with the gain factors. The results presented in table Q] correspond again to 
the most interesting case in which LISA has all its six laser links working. Notice 
that for all the parameters, the median values in the case of the FWF distributions 
show an improvement with respect to the RWF ones, which for some of them is more 
than one order of magnitude. The difference between the 10 % and 90 % levels of 
cumulative probability distributions provide the error intervals for 80 % of the source 
population. 

From table [U the reader can see the remarkable improvements for the 1O 7 M0 — 
10 7 M Q , and 10 7 M Q — 10 6 M Q cases in angular resolution and distance measurement: 
the angular resolution improves the median value by factors of 38 and 9.7, respectively; 
and the luminosity distance by 76 and 6.4, respectively. One should also notice that, in 
these two cases, those parameters were poorly determined using only the RWF. For the 
other sets of masses the averaged improvement in angular resolution and luminosity 
distance is more moderate, between 2 — 3 for Afijv, and 1.6 — 2.4 for ADl/Dl. In 
all cases the masses are determined much more accurate, even by several orders of 
magnitude in the case of /x, using the FWF. For the equal mass cases, the errors in 
M and \i are of the same order using the FWF. 

Another interesting way of looking at these probability distributions, is computing 
the percentage of fiducial sources that would have an error less than a certain number, 
and compare how these populations change when we use one or other waveform 
model. In table [5] we summarize these quantities for some of the most interesting 
parameters. In our analysis we pay special attention to the improvement of LISA 
angular resolution due to the inclusion of the higher harmonics in the waveform. One 
of LISA's objectives is to detect the signal during the inspiral phase alone and estimate 
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Table 2. Percentage of sources from Monte Carlo simulations of table [T] with a 
measurement error less than the given values. 
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its sky location accurately enough in order to issue warnings to the astronomical 
community of possible simultaneous observation of merger (and ringdown) with X- 
ray and optical observatories. The aim would be to identify the electromagnetic 
counterpart (host galaxy or galaxy cluster) and be able to disentangle the redshift 
parameter, since exiting science can be done by studying the luminosity distance 
redshift relation [III HI HI [H E] . 

Taking into account real field of view values of some of the existing or planned 
wide field electromagnetic instruments (see table 1 of [23j ) , we have studied the number 
of sources that could be observed by LISA with an angular resolution better than 
(2.5° x 2.5°), (1° x 1°) and (0.5° x 0.5°) corresponding to 1.9 x 10~ 3 , 3.0 x 10" 4 and 
7.6 x 10~ 5 srad, respectively. For the mass cases considered in table [H we have that 
using the FWF at least 20 % of the binary systems, located at z = 1 observed by 
LISA during the last year of coalescence, will have an angular resolution better than 
(1° x 1°). This represents 15 % more of sources than if we were using the RWF. 
Moreover, except for the very high mass systems M > 10 7 M@, we would have at 
least 14 % of sources resolved with an angular resolution better than 0.5° x 0.5° if 
we worked with FWF, while these numbers drop to less than 4 % for the RWF case. 
Therefore, we can conclude that the inclusion of the higher harmonics in the waveform 
is crucial in terms of being able or not to make observations of their electromagnetic 
counterpart. Concerning the other parameters in table [21 luminosity distance and 
masses, the errors are less than 10 % in almost all the cases using the FWF, while this 
is not the case with the RWF model. Notice that results depend a lot on the masses 
of the binary. 

In order to issue warnings and search for electromagnetic counterparts, we need 
LISA to have enough angular resolution some time before the merger occurs. For this 
reason, we are not only interested in LISA's angular resolution when we observe a 
source during the last year before merger, but also how the angular resolution evolves 
when we stop our observation some time before, and measure the importance of the 
FWF versus RWF. In figure [T] we study how the angular resolution is accumulated 
as a function of the lookback time, tjgco ~ tfi f° r different sources located at 
cos#jv = —0.383, 4>n — 2.82. Different sky locations have also been analyzed, see 
e.g. figure 14 in [THj for cos6>jv = —0.6, 4>n = 1- In this latter figure, we have that an 
angular resolution twice the final one is reached, respectively, 20 days, 20 hours and 
8 hours before merger, for binary systems of masses (10 b — 10 6 ) Mq, (10 7 — 10 6 ) M Q 
and (10 7 — 10 7 ) Mq. For the sky location considered here (figure [1]), LISA's angular 
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Figure 1. Time evolution of LISA'S angular resolution when we observe different 
supermassive black hole binary systems from one year before the merger up to 
tf. Lookback time is represented the axis of abscissas, so time evolves from right 
to left. Fiducial sources are located at redshift z = 1, with a common orientation 
defined by cos0^ = 0.2 and ifr^ = 3 and a sky location of cosSjv = —0.383 ; 
<f>N = 2.82. Solid lines correspond to FWF, and the dashed ones to RWF. 

resolution decreases more during the last days (hours) before merger, but we still 
are reaching almost the best angular resolution 12 hours before merger. The jumps 
in the evolution of LISA's angular resolution correspond to those times in which a 
new harmonic enters into the LISA band and it is related to the the low-frequency 
cut-off we have assumed for the instrument at 5 x 10 ~ 5 Hz. For equal masses, only 
the even multipoles of the orbital frequency contribute to the FWF, while for the 
unequal masses there is contribution from all the harmonics. For example, in the 
10 7 M© — 1O 7 M0 case the contribution of the 4th harmonic becomes relevant around 
10 days before coalescence (independently of the sky location considered) while the 
2nd harmonic rapidly increases the angular resolution 2 days before coalescence. For 
the unequal mass case 10 7 M Q - 10 6 M Q we clearly sec the contributions of the 2nd, 
3rd and 4th harmonics. 

Up to this point, we have considered LISA as two independent Michelson 
interferometers, using the long wavelength approximation, and we have provided the 
measurement errors for this two-detector case. But it is also interesting to study which 
is the impact on parameter estimation for a single interferometer. This is relevant in 
the hypothetical case that LISA would lose one of its arm links, in which case it would 
be interesting to know how LISA's parameter estimation gets reduced and whether 
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Table 3. Comparison between median values of the probability distributions 
when we see LISA as a single Michelson interferometer or as a combination of 
two independent ones. Gain factors are computed in the same way as in table [T] 
and by loss we mean the factor that the errors get multiplied by (and the SNR 
divided by) when we lose one of the LISA'S arms. 
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(c) mi = 
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the effect is the same when we use the RWF model or the FWF one. 

In table [3j we give the 50 % levels of the cumulative probability distributions 
working with both RWF and FWF, for the single detector and the two-detector case. 
As one would expect, each of the independent interferometers measures, in average, 
similar values of SNR, and since p tot = \J (p 1 ) 2 + (p 11 ) 2 , losing one of LISA's link 
means that the SNR gets reduced by a factor \f2 in all cases. For the measurement 
errors, we see that in most of the cases (especially for high mass systems) the loss 
factor (i.e. the factor that errors have to be multiplied by when LISA loses one of its 
links) is higher for the RWF case, which means that working with the FWF not only 
improves LISA's parameter estimation as we have seen before, but also reduces the 
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impact of losing one arm link. In any case, even working with the FWF model, we see 
that losing one arm, means reducing the angular resolution by more than one order of 
magnitude, thus we will lose any possibility of observing the potential electromagnetic 
counterpart of the merger. Losses in the other parameters are not so serious: errors 
in the distance get multiplied by a factor 2 — 5 and in masses by less than 2. 

4. Summary and outlook 

This work clearly shows that modeling the inspiral with the full post-Newtonian 
waveform, as compared to the restricted-PN one, improves in general the parameter 
estimation of supermassive black hole systems, due to the much greater richness 
of the waveform. One should notice that the results presented here are influenced 
by a number of assumptions associated to these observations, in particular, the 
instrumental and confusion noise models we have adopted, the post-Newtonian 
waveform model and how this is terminated, and the lower-frequency cut-off we 
have imposed. The parameter estimation code we use has been validated by the 
LISA science performance evaluation taskforce [24| . This new taskforce, created at 
the September'07 LISA International Science Team meeting, is currently extending 
these results to include spin-induced precession into the full waveform and compare 
LISA's science reach for different mission configurations, based also on astrophysically 
motivated source distributions [25j . It will be very interesting to extend these analyses 
using different post-Newtonian models and also to study parameter estimation for all 
three stages of the signal, including the merger and the ringdown, and use currently 
available phcnomenological and numerical waveforms for that. 
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